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Global  Measurements  of 
Stratospheric  Mountain  Waves 
from  Space 

Stephen  D.  Eckermann1  and  Peter  Preusse2 

Temperatures  acquired  by  the  Cryogenic  Infrared  Spectrometers  and  Tele¬ 
scopes  for  the  Atmosphere  (CRISTA)  during  shuttle  mission  STS-66  have  pro¬ 
vided  measurements  of  stratospheric  mountain  waves  from  space.  Large-am¬ 
plitude,  long-wavelength  mountain  waves  at  heights  of  15  to  30  kilometers 
above  the  southern  Andes  Mountains  were  observed  and  characterized,  with 
vigorous  wave  breaking  inferred  above  30  kilometers.  Mountain  waves  also 
occurred  throughout  the  stratosphere  (15  to  45  kilometers)  over  a  broad 
mountainous  region  of  central  Eurasia.  The  global  distribution  of  mountain 
wave  activity  accords  well  with  predictions  from  a  mountain  wave  model.  The 
findings  demonstrate  that  satellites  can  provide  the  global  data  needed  to 
improve  mountain  wave  parameterizations  and  hence  global  climate  and  fore¬ 
cast  models. 


Gravity  waves  in  the  lower  atmosphere, 
forced  by  flow  over  mountains,  have  been 
observed  and  modeled  for  many  years  (1). 
Mountain  waves  can  also  propagate  into  the 
stratosphere.  Drag  produced  by  breaking 
stratospheric  mountain  waves  is  essential  for 
accurate  climate  and  forecast  models  (2).  The 
turbulence  generated  by  stratospheric  moun¬ 
tain  wave  breaking  poses  a  major  safety  haz¬ 
ard  to  high-altitude  aircraft  (3-5)  and  influ¬ 
ences  stratosphere-troposphere  exchange  of 
important  trace  gases  (6).  Mountain  wave 
ascent  triggers  formation  of  polar  stratospher¬ 
ic  clouds  (PSCs)  (7).  It  has  been  argued  that 
mountain  wave  PSCs  reduce  Arctic  ozone 
concentrations  by  activating  chlorine  in 
stratospheric  air  that  flows  through  them  (8). 

Computational  limits  have  prevented  global 
troposphere-stratosphere  models  from  operat¬ 
ing  at  a  spatial  resolution  fine  enough  to  gen¬ 
erate  mountain  waves  in  the  stratosphere,  so 
their  effects  must  be  parameterized  (2).  The 
spatial  resolution  of  stratospheric  sounding  sat¬ 
ellites  has  generally  been  too  coarse  to  resolve 
mountain  waves.  Detailed  measurements  have 
been  limited  largely  to  scattered  ground-based 
soundings  (9),  some  regional  campaigns  (3, 
10),  and  chance  intercepts  by  high-altitude  air¬ 
craft  (4,  11).  Thus,  global  data  are  lacking. 

Recently,  a  few  high-resolution  strato¬ 
spheric  satellite  instruments  have  detected 
gravity  waves  (12).  However,  instrumental 
effects  have  complicated  interpretation  of 
these  gravity-wave  signals  (13).  In  particular, 
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the  data  have  shown  little  clear  evidence  of 
mountain  waves,  despite  suggestions  that 
mountains  should  be  a  strong  source  of 
stratospheric  gravity  waves  (2-4,  9-11). 

On  4  November  1994,  the  CRISTA-SPAS 
(Shuttle  Pallet  Satellite)  experiment  was  de¬ 
ployed  into  orbit  by  the  space  shuttle  Atlantis 
(Fig.  1A).  CRISTA  measured  infrared  emission 
spectra  (4  to  71  pm)  from  atmospheric  limb 
scans  with  the  use  of  three  telescopes  that 
viewed  closely  separated  volumes  of  air  (14). 
From  the  CO-,  Q-branch  emission  at  — 12.6  pm, 
vertical  temperature  profiles  T(z)  with  a  preci¬ 
sion  of  ~0.5  K  were  derived  every  1.5  km  in 
altitude,  z,  throughout  the  stratosphere  (15,  16). 
Although  short  spatial  scales  along  the  line  of 
sight  are  smeared  out,  medium  scale  structure  is 
resolved.  Theoretical  radiance  calculations 
show  that  the  temperature  data  have  >50% 
visibility  to  gravity  waves  with  wavelengths  S3 
to  5  km  vertically  and  a  100  to  200  km  along 
the  limb  (17).  Each  telescope  provided  a  profile 
every  200  or  400  km  along  the  orbital  track, 
depending  on  the  measurement  mode  (18),  for  a 
week. 

Small-scale  temperature  fluctuations  T'(z) 
were  isolated  from  global  temperature  data  with 
a  wavenumber  0  to  6  Kalman  filter.  Vertical 
wavelengths  and  amplitudes  T  of  dominant 
oscillations  in  each  T'(z)  profile  were  deter¬ 
mined  with  spectral  peak  detection  and  sliding 
harmonic  fits,  respectively.  These  parameters 
are  shown  (Fig.  IB)  for  three  successive  T'(z) 
profiles  (labeled  1,  2,  and  3)  acquired  over  the 
southern  Andes  Mountains  (Fig.  1C)  on  6  No¬ 
vember  1994.  All  three  profiles  show  similar¬ 
looking  large-amplitude  wave  oscillations  at  15 
to  30  km.  The  photo  in  Fig.  1  A,  taken  over  this 
region  2  days  earlier,  shows  banded  mountain 
wave  clouds  downstream  of  the  Andes.  This 
suggests  that  the  oscillations  in  Fig.  IB  may  be 
stratospheric  mountain  waves. 


A  long,  stationary  plane  mountain  wave 
emanating  from  the  Andes  should  have  a 
vertical  wavelength  given  by  the  hydrostatic 
relation 

__  2nU 

theory  (1) 

where  U  is  the  local  horizontal  wind  speed 
across  the  Andean  Ridge  and  N  is  the  buoyancy 
frequency  (1).  At  15  to  30  km,  U  =  20  to  23  m 
s_1  (Fig.  ID)  and  N  =  0.020  rad  s  ',  which 
yields  (A..),^  ^  6  to  7  km.  The  observed 

vertical  wavelength  of  the  measured  wave 
oscillations  in  Fig.  IB  is  ~6.5  km  at  15  to  30 
km.  The  close  agreement  between  observed  and 
theoretical  vertical  wavelengths  is  strong  evi¬ 
dence  that  these  measured  wave  oscillations  are 
stratospheric  mountain  waves.  Successive  pro¬ 
files  are  separated  horizontally  by  200  km  (Fig. 
1C)  and  phase  shifted  by  180°  (Fig.  IB).  This 
implies  a  horizontal  wavelength  \h  =  400/;?  km 
for  this  wave,  where  n  is  an  odd  integer.  Given 
the  sensitivity  of  CRISTA  to  \h  >  100  to  200 
km,  then  \h  must  be  either  130  or  400  km. 

Nonbreaking  mountain  waves  grow  in  am¬ 
plitude  with  height  approximately  as  exp (z/ 
2Hp).  The  density  scale  height  H  ~  1  km,  and 
so  T  should  increase  by  a  factor  of  3  over  the 
15-  to  30-km  height  range,  as  observed  in  Fig. 
IB.  Thus,  the  mountain  wave  propagates  with¬ 
out  breaking  from  15  to  30  km.  Above  30  km, 
the  oscillations  in  Fig.  IB  abruptly  attenuate 
and  disappear.  There  are  two  plausible  reasons 
for  this.  Because  winds  weaken  rapidly  above 
this  level  (Fig.  ID),  decreases  according  to 
Eq.  1,  and  the  wave  soon  becomes  too  short  for 
CRISTA  to  resolve.  Additionally,  however,  the 
wave  is  likely  to  break  above  30  km.  Mountain 
waves  advect  heavy  air  over  lighter  air  and 
overturn  when  amplitudes  exceed  the  breaking 
threshold 


where  is  the  adiabatic  lapse  rate  (9.8  K 
km~ ')  and  T  is  background  temperature  (19).  At 
15  to  30  km,  the  mountain  wave-breaking  am¬ 
plitude  7’brcak  10  K.  This  value  is  consistent 
with  the  observed  nondissipative  growth  in  f 
from  15  to  30  km,  because  f  <  tbteak  at  these 
heights.  At  30  km,  however,  t  =  7  K  (Fig.  IB) 
and  is  probably  larger  in  reality,  given  some 
observational  smearing.  Indeed,  retrieval  calcu¬ 
lations  indicate  that,  for  a  wave  of  X.  =  6.5  km 
and  \h  =  400  km,  the  amplitude  T  measured  by 
CRISTA  is  only  ~60%  of  its  true  value  (17). 
Thus,  this  wave  is  at  or  near  its  breaking  ampli¬ 
tude  at  30  km.  Above  30  km,  Eq.  2  shows  that 
rbreak  decreases  because  of  the  weakening  winds 
in  Fig.  ID,  making  wave  breaking  and  ampli¬ 
tude  attenuation  even  more  probable  and  severe. 

The  influence  of  a  breaking  mountain  wave 
on  the  atmosphere  is  determined  by  its  vertical 
flux  of  horizontal  momentum  density,  the  so- 
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called  Eliassen-Palm  (EP)  flux,  which  can  be 
computed  from  measured  wave  parameters 
(20).  With  the  use  of  Eq.  2  and  linear  saturation 
theory  (2,  19),  wave  breaking  above  30  km 
yields  large  EP  flux  divergences  that  imply 
mean-flow  accelerations  of  up  to  —(50  to  200) 
m  s-1  day-1.  These  decelerative  forces  on  the 
eastward  stratospheric  flow  in  Fig.  ID  are  one 
to  two  orders  of  magnitude  larger  than  zonal- 
mean  climatological  estimates  of  small-scale 
wave  forcing  at  40°  to  50°S  (21).  Such  large 
values  imply  vigorous  local  wave  breaking, 
drag,  turbulence  production,  and  mixing  at  30 
to  40  km.  The  wave  is  likely  to  dissipate  en¬ 
tirely  below  43  km,  where  a  critical  level  (U  = 
0:  Fig.  ID)  causes  both  and  T  to  vanish  (Eqs. 
1  and  2)  and  prevents  further  penetration  of 
mountain  wave  energy  ( 1 ). 

How  common  are  large  mountain  wave 
events?  Figure  2  compares  and  (k-)theory  at 
z  =  25  km  based  on  T'(z)  data  for  all  seven 
mission  days  at  40°  to  50°S  over  South  Amer¬ 
ica.  Again,  there  is  very  good  correlation  be¬ 
tween  theoretical  and  measured  vertical  wave- 
lengths.  Amplitudes  measured  on  four  succes¬ 
sive  days  are  plotted  in  Fig.  3,  A  to  D.  Under¬ 
lying  plots  show  corresponding  predictions, 
with  the  Naval  Research  Laboratory  Mountain 
Wave  Forecast  Model  (MWFM)  ( 4 ,  22)  inter¬ 
faced  to  assimilated  regional  winds  and  temper¬ 
atures  from  NASA’s  Data  Assimilation  Office 
(DAO)  (Fig.  3,  E  to  G)  (23).  The  MWFM 
results  capture  much  of  the  day-to-day  variabil¬ 
ity.  Figures  2  and  3  reveal  that  mountain  waves 
were  present  in  this  region  of  the  stratosphere 
throughout  the  mission. 

These  long  plane  stratospheric  mountain 
waves  resemble  those  observed  by  aircraft  over 
the  Scandinavian  Ridge  (8,  24).  Rapid  cooling 
of  air  by  the  latter  waves  produces  PSCs,  which 
in  turn  destroy  Arctic  ozone  (8).  Although  the 
Andean  mountain  waves  observed  here  also 
produce  large  temperature  decreases  (Fig.  3), 
they  are  probably  too  far  north  to  yield  temper¬ 
atures  cold  enough  to  form  PSCs.  Nevertheless, 
as  the  edge  of  the  Antarctic  vortex  moves  back 
and  forth  across  the  tip  of  South  America  (25), 
strong  turbulence  generated  by  these  mountain 
waves  will  mix  tongues  of  ozone-depleted  vor¬ 
tex  air  with  ozone-rich  air  from  mid-latitudes, 
thus  diluting  local  ozone  concentrations  overall. 
The  large  EP  flux  divergences  could  also  mod¬ 
ify  the  subsequent  evolution  of  the  vortex  (26). 

Are  stratospheric  mountain  waves  evident 
elsewhere  in  the  CRISTA  data?  T  amplitudes 
over  the  Northern  Hemisphere  at  25,  35,  and  45 
km  on  9  November  1994  are  plotted  in  Fig.  4, 
A  to  C.  Mountain  wave  amplitudes  T  predicted 
by  MWFM  on  this  day  are  shown  in  Fig.  4,  D 
and  F.  Although  the  peak  amplitudes  differ 
(27),  observations  and  model  results  show  very 
similar  geographical  distributions.  Most  notable 
is  a  broad  region  of  enhanced  amplitudes  over 
central  Eurasia.  The  activity  is  produced  by 
mountain  waves  from  the  Altai  Mountains, 


Fig.  1.  (A)  Astronaut  photograph  over  southern  South  America  showing  deployment  of  CRISTA- 
SPAS  from  Atlantis,  4  November  1994,  at  12:19  UT.  The  cylindrical  assembly  on  the  satellite 
houses  CRISTA.  The  east  coast  of  Argentina  and  the  Atlantic  Ocean  are  visible  to  the  right;  banded 
mountain  wave  clouds  downstream  of  the  Andes  are  visible  to  the  lower  left.  (B)  Temperature 
perturbations  T'(z)  from  three  successive  CRISTA  scans  (labeled  1,  2,  and  3)  on  6  November  1994, 
at  —06:00  UT.  (C)  Topographic  elevations  for  southern  South  America  (linear  color  scale).  Squares 
labeled  1,  2,  and  3  show  locations  of  the  profiles  in  (B).  The  nadir  point  and  orbital  motion  of  the 
shuttle  in  (A)  are  shown  with  a  black  and  white  circle  and  black  arrow,  respectively.  (D)  Horizontal 
winds  orthogonal  to  the  Andean  Ridge,  from  DAO  assimilated  data  (23)  at  42°S,  69°W  on  6 
November  1994,  06:00  UT. 


Fig.  2.  Predicted  moun¬ 
tain  wave  vertical  wave¬ 
lengths  (\z)theo^  versus 
measured  wavelengths 
\z  for  all  profiles  over 
South  America  between 
40°  and  50°S  in  the  20- 
to  30-km  height  range. 

In  the  calculations  of 
(^z)theory'  U  was  approx¬ 
imated  by  the  average 
total  wind  speed,  be¬ 
cause  stratospheric  flow 
was  generally  orthogo¬ 
nal  to  the  Andean  Ridge 
over  the  mission  period. 

Wind  speeds  were  taken 
from  United  Kingdom 
Meteorological  Office 
assimilated  data  (23)  for 
the  days  and  times  in 
question,  with  gridded 
values  interpolated  to 

the  specific  profile  location  and  averaged  from  20  to  30  km.  Regions  of  largest  amplitude-weighted  point 
density  are  shaded. 
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Fig.  3.  f  amplitudes 
over  southern  South 
America  from  6  to  9 
November  1994  (A  to 
D,  respectively),  plot¬ 
ted  as  color  pixels  at  the 
profile  location.  Raw 
values  were  scaled  ac¬ 
cording  to  their  verti¬ 
cal  wavelengths  to 
correct  for  observa¬ 
tional  smearing  (17). 
Beneath  each  panel  is 
a  corresponding  MWFM 
prediction  based  on 
DAO  assimilated  winds 
and  temperatures  at 
12:00  UT  for  these 
dates  (E  to  H).  Modeled 
waves  with  \z  <  5  km 
were  not  plotted  be¬ 
cause  CRISTA  is  not 
sensitive  to  them.  The 
color  scales  are  linear. 
Pixels  are  plotted  in  or¬ 
der  of  ascending  ampli¬ 
tude  and  are  sized  to 
maximize  coverage  and 
visibility  while  minimiz¬ 
ing  pixel  overlap. 


T 
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Fig.  4.  (A  to  C)  Polar 
orthographic  maps  of 
amplitudes  f  from 
CRISTA  profiles  over  the 
Northern  Hemisphere 
on  9  November  1 994  at 
25,  35,  and  45  km,  re¬ 
spectively.  Values  were 
corrected  for  observa¬ 
tional  smearing  as  in 
Fig.  3.  (D  and  F)  Corre¬ 
sponding  MWFM  am¬ 
plitudes  f  at  25  and  45 
km,  respectively,  based 
on  DAO  assimilated 
winds  and  tempera¬ 
tures  at  12:00  UT  on 
this  date.  As  in  Fig.  3, 
modeled  waves  with  \z 
<  5  km  were  omitted. 
(E)  Topographic  eleva¬ 
tion  in  the  Northern 
Hemisphere.  All  color 
scales  are  linear.  Pixels 
are  plotted  in  order  of 
ascending  amplitude. 


CRISTA  T :  25  km  CRISTA  T :  35  km  CRISTA  T :  45  km 
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Sayan  Ranges,  and  Tian  Shan,  which  skirt  the 
border  regions  of  Russia,  Mongolia,  China,  and 
Kazakhstan  (Fig.  4E).  These  mountain  waves 
are  present  at  all  heights,  reaching  the  upper¬ 
most  levels  of  the  stratosphere  (Fig.  4C).  There 
is  also  evidence  of  more  isolated  mountain 
wave  activity  over  the  Alps  and  the  Alaska- 
Yukon  region.  Mountain  waves  are  largely 
absent  in  both  the  observations  and  model 
results  over  other  mountainous  regions, 
such  as  western  North  America,  the  Hima¬ 
layas,  Greenland,  and  Scandinavia.  This 
absence  occurs  in  the  MWFM  results  be¬ 
cause  regional  winds  inhibit  mountain 
wave  propagation  into  the  stratosphere. 

Despite  the  basic  similarities,  specific  dif¬ 
ferences  between  the  CRISTA  data  and 
MWFM  results  are  also  evident  in  Figs.  3  and 
4.  Shortcomings  in  the  model  may  be  impli¬ 
cated  here,  highlighting  the  potential  value  of 
observations  such  as  these  in  improving  glob¬ 
al  mountain  wave  parameterizations. 

Figure  4C  also  shows  large  t  amplitudes  at 
z  —  45  km  at  the  northernmost  latitudes,  which 
cannot  be  explained  in  terms  of  polar  mountain 
waves  (Fig.  4F).  This  observation  is  consistent 
with  large  gravity-wave  amplitudes  observed 
previously  in  ground-based  soundings  of  the 
high-latitude  winter  upper  stratosphere  (28). 
The  feature  has  been  attributed  to  enhanced 
exp(z/2Hp)  amplitude  growth  of  weaker  “back¬ 
ground”  gravity-wave  activity,  due  to  the  tem¬ 
perature  structure  of  the  winter  polar  strato¬ 
sphere  (13,  29),  although  longitudinal  asymme¬ 
tries  in  stratospheric  winds  can  modulate  the 
effect  substantially  (30).  This  indicates  that 
CRISTA  measured  gravity  waves  from  sources 
other  than  mountains  as  well. 
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tribution  tends  to  be  controlled  largely  by  cli¬ 
mate  (8,  9),  and  surface  property  changes  can 
affect  climate  by  modifying  the  atmospheric 
energy  and  water  budget  (10-13),  it  is  reason¬ 
able  to  propose  that  dynamic  vegetation-climate 
interaction  might  influence  decadal  climate 
variability  substantially  in  a  climatically  sensi¬ 
tive  zone  such  as  the  Sahel.  We  tested  this 


Enhancement  of  Interdecadal 
Climate  Variability  in  the  Sahel 
by  Vegetation  Interaction 

Ning  Zeng,1*  J.  David  Neelin,1  K.-M.  Lau,2  Compton  J.  Tucker2 

The  role  of  naturally  varying  vegetation  in  influencing  the  climate  variability 
in  the  West  African  Sahel  is  explored  in  a  coupled  atmosphere-land-vegetation 
model.  The  Sahel  rainfall  variability  is  influenced  by  sea-surface  temperature 
variations  in  the  oceans.  Land-surface  feedback  is  found  to  increase  this  vari¬ 
ability  both  on  interannual  and  interdecadal  time  scales.  Interactive  vegetation 
enhances  the  interdecadal  variation  substantially  but  can  reduce  year-to-year 
variability  because  of  a  phase  lag  introduced  by  the  relatively  slow  vegetation 
adjustment  time.  Variations  in  vegetation  accompany  the  changes  in  rainfall, 
in  particular  the  multidecadal  drying  trend  from  the  1950s  to  the  1980s. 
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